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ABSTRACT 

In this paper the elasticity problan for a laminated thick plate 
which consists of twc bonded dissimilar layers and which contains a 
circular hole Is considered. The problem Is formulated for arbitrary 
axlsynmetrlc tractions on the hole surface by using the Love strain 
function. Through the expansion of the boundary conditions Into Fourier 
series the problem Is reduced to an Infinite system of algebraic equa- 
tions which Is solved by the method of reduction. Of particular Inter- 
est In the problem are the stresses along the Interface as they relate 
to the question of delamination failure of the composite plate. These 
stresses are calculated and are observed to become unbounded at the hole 
boundary. An approximate treatment of the singular behavior of the stress 
state Is presented and the stress Intensity factors are calculated. 


1. INTRODUCTION 

In this paper the elasticity problem for a thick plate which con- 
sists of two bonded dissimilar homogeneous layers Is considered. It is 
assumed that the plate Is Infinite, contains a circular hole, and Is 
subjected to axl symmetric external loads. Even though the problem as 
stated may have some applications, from a practical viewpoint the Im- 
portant problem Is that of a laminated plate containing a circular hole 
and subjected to uniaxial membrane or bending loads away from the hole 
region. The latter problem has Important applications In the analysis 
of delamination failure of perforated multilayered plate and shell 
structures. In such structures the Interface stresses are known to 
have a power singularity which greatly enhances the possibility of 

*Th1s work was supported by NSF under the Grant CM E-7809737 and by NASA- 
Langley under the Grant NGR-39-007-011 

Department of Mechanical Engineering and Mechanics, Drexel University, 
Philadelphia, Pa. 
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delamlnatlon failure [1.2,3]. In a laminated plate under general load- 
ing conditions one may always separate a homogeneous solution and reduce 
the problem to a perturbation problem In which the self-equilibrating 
tractions on the hole surface are the only external loads. By expanding 
these tractions Into Fourier series In e one may further separate the 
problem Into Its simpler components In the Independent variables r and 
z only. Thus, the axl symmetric problem treated In this paper may also 
be considered as the first component of the general three-dimensional 
plate problem. 

The three-dimensional elasticity problem for laminated plates con- 
taining a hole does not seem to hav« been considered before (see, for 
example, [4] for a recent review). The existing solutions are mostly 
based on numerical techniques and are generally highly approximate [5-7]. 
The circular hole problem for a homogeneous thick plate was considered 
In [8] and [9]. The technique developed in [9] will be used to solve the 
laminated plate problem considered In this paper. The related problem 
of a layered semi-infinite medium (I.e., the limiting case of the hole 
problem In which the hole radius Is infinite) was considered in [10] 
where the method of singular Integral equations was used to solve the 
problem. 

2. THE HOMOGENEOUS SOLUTION 

In this section we give the results of some elementary solutions 
for a laminated plate without a hole which Is subjected to certain uni- 
form loading conditions at Infinity. First we assume that the composite 
plate which consists of two layers with the elastic constants E-j, and 
Eg, vg, and the thicknesses h^ and hg Is subjected to an average radial 
membrane stress at r*« and Is constrained to remain flat upon defor- 
mations. Thus, defining 

o^^l(r,e,z) = a■^ , a^^g(r,0.z) * Og , (1) 

we have 
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°eel * 


f 




( 2 ) 


and all the remaining stress components are zero. From 


c,h, + Cjhj . o„(h,+Kj) . • £^^2 


( 3 ) 


It then follows that 

(l-v2)m(l+k) 

'^1 ’ “"o T-v,+OT:^1 » 

(l-vi)(l+k) 

“ % r-v^+ikTr-TgT ’ 

where 

m * E-j/Eg , k * h-j/hg . 

Next we consider the composite plate which Is uniformly loaded by 
an average stress In x-direction only. 

Defining now 

Oxxi(^»y*^) * ~ ^2 * 

and again assuming that the plate Is constrained to remain flat, the 
nonzero stress components op ag. determined from 

°lhl ®2^2 “ '’o^h'^^2^ * Sxl “ Sx2 » 

"yylh * °yyZ^2 ' ° * Vyl * Vy2 * 

Solving (7) we have 


(4a.b) 


(5) 
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ofQ(l+k)[m(l-v^v2) + m2k(l-v2)] 

^ 1 -v| + 2mk(l-v.j V2) 

02 " <JQ(l+k) - o^k , 

a, - mo2 

®yy1 * V.J + mkv2 * *^yy2 * ’*^‘^yyl ‘ 

Referring to the cylindrical coordinates the stress states in the 
layers 1 and 2 may be expressed as 




0 .+0 


u • G * *0f * 

= cos 20 , 


eej 


Vej 


- J sin 2e , 


°22j ” °r2j ■ “62j " ° 


(9»-f) 


Solutions similar to that given by (4) and (8) may also be obtained 
for other types of uniform external loads such as bending and thermally 
induced loading. 


3. SOLUTION OF THE AXISYMMETRIC PROBLEM 

The problem of a laminated plate containing a circular hole with a 
radius a may now be solved by superimposing on the homogeneous solutions 
found in the previous section a perturbation solution in which the trac- 
tions -o^^j(a,0,z), (a»r,0,2; j=l,2; O<0<2 it, -h2<z<h.|) acting on the hole 
boundary are the only external loads, where the stress state 
(e,a “ r,0,z; j®l,2) is given by the homogeneous solution. From, for 
example. (1), (4), (8), and (9) it may be seen that the simplest such 
perturbation problem is an axi symmetric problem corresponding to the 
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axlsymnetrlcally loaded plate or to the first part of the unldlrectlon- 
ally loaded plate (I.e., to the e>1ndependent part of the solution given 
hy (9)}. Thus, In this section Me will consider the axl symmetric prob« 
lem for the composite plate subjected to the following boundary condi- 
tions (Figure 1): 

Orrj(a»2) • Oj(z) . (10a,b) 


where and tj are known tractions acting on the layers 1 and 2, re- 
spectively. Note that the solution Is Independent of e and the shear 
stressBOggj = 0, (8*r,z; j=l,2) everywhere. 

To formulate the problem the technique described In [9] will be 
used. Because of axlsymmetry. It Is sufficient to use the z-component 
of the Galerkln vector only which is nothing but the Love strain func- 
tion Z(r,z) [11,12]. In addition to the surface tractions given by 
(10) the composite plate Is subjected to the following homogeneous 
boundary, continuity, and regularity conditions (Figure 1): 


<yrzi(>^»f’l) “ 0 » 

‘^rz2^’**'^2^ “ ° * <^zz2^^*"^2^ * ° * (a<r<«) , (12) 

^rzl^^*^^ * * *^zzl^***^^ * ®2 z 2^***^^ * (a<r<*) ,(13) 

Upl(i*.0) ’ “ 022^''*®^ ’ (a<r<») , (14) 

”rrj(-, 2 ) - 0 . 

Let and be the Love strain functions for layers 1 and 2, re- 
spectively (Figure 1). In the absence of body forces Z-| and Zg satisfy 

72v2Zj(r,2) • 0 . + + . (j-1.2) . (16) 
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and the displacements and stresses are given by 


2vur • - . 2uu^ - [2(l-v)v2 - ^]Z ; (17) 

°rr “ 32 3r^^^ * ‘^ee * 3z r ?r^^ * 

®Z2 ’ « C(2-v)v2 - ^]Z . ^rz “ ^ [(l-v)v2 - ^]z . (18) 


Looking for a series type solution In z and taking Into account the 
regularity conditions (15), the solution of (16) may be expressed as 

n*i 

Substituting from (17)-(20) Into the boundary and continuity con- 
ditions (11)-(14) and observing that these conditions must be satisfied 
for all values of r In a<r<®, after some manipulations. It may be shown 
that part of the resulting algebraic system which conwalns only the 
unknowns and is separated and Is overdetermined. 

Hence for the conditions (11)-(14) to be satisfied it Is necessary that 


'3n °3n “ ''3n "'3n 


(21) 


Considering (21) from (17)-(20) the displacements and stresses In the 
composite plate may now be expressed as follows: 
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“n * ^7 ‘ “n''l<“n'■>C^n * V“n * 

m 

■ nil * ®5n*“n ' *Sn3*'%*> • 

“^' ■ ^‘nil '<"“n«o<“n'‘>C(-2*^''l>®5n * *1n 

* *5n^n3 * j, “S “s»n*'<o(“n'-)C<2-^'’l)*5n 

*®ln*V“n]> • 

“rrl ' J, “n<®lt.[''o(“n'') '^<V>3 

* ®5n^C“n''o<“n''> * 7 l<l(an'')] 


* "l<“n''> * n+2v,)K„(c.„r)])s1na„2 
+ i|C,(a„r)]*B5„[(H2v,)K„(,„r) 

*^"l<V>l> “*V • 

"eel ■ - J, “J C^ln * ®5n*“n * *5nl 

* 2'’l“n*5n’^o(“n''>> * J, “S “*“n« “r" C*ln 

* *5n^°n *■ ®5n^ ' ^''l“n®5n''o*“n'''* > 


"22I • - J, “n "o<“n'^)[®ln * ®5n"“n * (1-2v,)A5„]slna„2 


s 
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00 

"rzl ■ - „I, ♦ A,„ * 

00 

“rt ' ^ < J, l'l(»„'-)[C,„ * C5„z«„ ♦ D5„]cosa„A 

oo 

■ „l, “n * V“n ' • 

“z2 <‘1 s1na„zK„(a„r)[(-2+4v2)05„ 

00 

* '=ln * <:5n^“n3 * J, “S«>S‘‘n^Ko(“n'')[(2-4v2)C5„ 

* “in * “Sn4“n^> ’ 

00 

V2 * J, '‘n'“lnt''o<“n'') * 

* “5n4f“n''o('‘n'‘) ^ ? I<1 (“„'•)] 


(26) 


(27) 


(2B) 


(29) 


■ “sn^^ '^l(“n''* * ('*7'’2)Ko(»„'')]tsina„z 

00 

■ „^1 “n<“lnC^(»„'') * ^ '<,(»„'•)] ♦ C5„z[.„K„(a„r) 

* r "l*“n''*^ * “snf*’*^''2*''o*“n'"' * ^ ''l*“n'^*3 )cosa„z. (30) 


’992 ■ - “5 s1"«„a(-S: 2- C“ln * “sn^-, 
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• “n "”n*< 

C'^ln * * “Sn^ - V.“srVV>> • <3D 

0» 

”uZ ■ • J, “n ^)<“n''>C»l„ * V“r. * ('-2''2)C5„]s1n<.„z 

* „I, “n "o<“n''>C'ln * '=5n“n * (^vj-DDj^Jeosa^z . (32) 

OQ 

■ ■ „^i “n'^l'“n''>*'%^t2''2'>5n * ^In * 'sn^n^ 

00 

■ 1 °^l<“n''>‘=““n*[-2''2C5n * “in * “sn^’n^ • (33) 


Substituting from (22)-(33) Into the boundary and continuity conditions 
(n)-(14) we obtain the following system of homogeneous algebraic equa- 
tions: 


-A^nSinanh^ + [2viC0Sa„h^ - «nh^s1na„h,]A5„ 

-Bincosanh^ - [2v,s1no„h^ + On^'iCOSa^h^lBg^ « 0 , (34) 

A^nCOSanh^ + [(2 vt 1 )s1na^h^ + ‘ ®ln®^%^ 

+ C(2vyl)cosa^h^ - “nh^sIno^h^lBgj, = 0 . (35) 

CinSln«^h2 + [2v2COs«„h2 - «n^2S^%^2^^5n 
- D^„C0Sa^h2 + [2v2Slna^h2 + “nh2COSa^h2]Dg^ » 0 , (36) 

^■]n‘^°®“n^2 ^CO-2v2)s1na^h2 - «nh2COSa^h2]C5^ 
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■*’ * C(2v2-1)coso^h2 - ° * 

- + 2v2(^'D,n ' ° * 

^ * ^In * • 0 . 

®^n * ®®5n ■ ^In ' ‘^Sn * ° ’ 

6(2-4v^)A 5^ + - ^2-4v2)Cg^ " hn ’ ^ * 

where 


(37) 

(38) 

(39) 
(^ 0 ) 
(41) 


S * . (42) 

To obtain a non-trivial solution for the system of algebraic equations 
(34)-(42), the determinant of coefficients must vanish, giving 


4(a„) « 0 (43) 

The characteristic equation (43) gives the eigenvalues (n»l,2,...) 
of the problem. A close examination of the roots of (43) shows that If 
Is a root so are -Op, Op, and -Op. Therefore In solving the problem It Is 
sufficient to use the roots In the first quadrant and consider the real 
part of the solution only. Further examination of the roots Indicates that 
for dissimilar materials generally there Is only one positive real root and 
the remaining roots are all complex. Furthermore, the complex roots form 
two distinct sequences In the first quadrant which greatly facilitates 
their numerical evaluation. When the elastic constants of the layers 1 
and 2 become equal (I.e., for a homogeneous layer), the real root disap- 
pears and the two sequences of complex eigenvalues become the roots of 
the following characteristic equations: 
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sin2x^ * \ * 0 , sin2x^ -2Xj, -0 . 

Equations (44) are known to correspond to the extension and bending prob- 
lems for a homogeneous thick plate containing a circular hole [6*9]. 

It Is found that a^*0 Is also a root of (43). Therefore, a partic- 
ular solution must be added to that given by (19) and (20) to account for 
the zero eigenvalue. Considered as sir.jle "plates" the layers are sub- 
jected to stretching, bending, and transverse shear. The particular 
solutions Z^, and must therefore exhibit the characteristics of all 
three modes of loading. Thus 

Z® ■ N^z*n r + M^z‘tn r + P^r2jrr+ Q.j(z2-rV2) , 

Z 2 • NgZ £n r + M 2 Z^£nr+ Pgr^sn r . (45a, b) 


where the first terms In each expression correspond to stretching, and 
the next two terms to combined bending and transverse shear In the Indi- 
vidual plates. The term Q^(z2-rV2) corresponds to a rigid body trans- 
lation In z-directlon and Is added to (45a) to Insure continuity of 
displacements at the Interface. The constants M<|, Ni, P^, N 2 , and 

P 2 which appear In (45) are not Independent. By using expressions (17) 
and (18) which relate the displacements and the stresses to the Love 
strain function, all field quantities can be written In terms of these 
constants. Then, by applying the boundary and continuity conditions 
(11)-(14), after some lengthy algebra and after redefining the constants 
we obtain: 


’''I 

Z° » DpZtn r ♦ r + r + Q^(z2-r2/2) 


^ ■ BO^jtn r + bAj,z*ui r ♦ 6 r 


(46a,b) 


Note that Z? and Z? are of the form f(r) g(r)h(z) + m(z) suggested 
in [93. 2 
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The stress field generated by (46) can then be expressed as: 

Vl ■ ^ ^ • 

° * “^rzl “ ^ ’ (47a-d) 

a°r2 “ * *^902 ' ' * 

4zZ ' ° ’ 4z2 ' ° 

In the perturbation problen the stress states in layers 1 and 2 are 
obtained by adding the respective stress components given by (24)-(27), 
(30)-(33), (47) and (48). Thus, the problem is reduced to one of determ- 
ining the unknown constants A^, D^, Ag^, Bg^, Cg^, 
and Dg^, (n*l,2,...). However, from (34)-(42) it is clear that the homo- 
geneous system contains only one arbitrary constant for each eigenvalue 
«j,, (n=l,2,...). For example, one may assume that A^^, (n»l,2..) is the 
only unknown in (34)-(42) and the remaining seven unknowns Ag^»-**»D 5 ^ 
may be expressed in terms of A.j^ after solving the related eigenvalue 
problem. The unknown constants A^, D^, and A^^^, (n=l,2,..) are then 
determined from the boundary conditions (10). To do this, we first 
substitute from the expressions (24), (47a), (30), (48a) and (27), (47d), 
(33), (48d) into (10a) and (10b), respectively. In the resulting equa- 
tions by expanding both sides into a series of an appropriate system of 
orthogonal functions in -h 2 <z<h-| and by matching the coefficients we 
obtain a linear system of algebraic equations to determine the unknown 
coefficients A^, D^, and A^^. The algebraic system is infinite and may 
be solved by the method of reduction. 

If we use the first N+1 functions of a real orthogonal system, the 
(real parts of) conditions (10) would give 2N+2 equations. On the other 
hand, since A^, 0^, A^^ (corresponding to the real eigenvalue a-|) are 
real and complex, truncating the series (24), (27), (30) 
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and (33) at the Nth term we would have 2N+1 real unknowns. However, 

It can be shown that this discrepancy disappears If one selects an 
orthogonal system In which the first function Is a constant. Thus, 

If we substitute from (27), (47d), (33), and (48d) Into 

f'o^dz . f » jdz + r'Vzl"" • <'•’> 

ihg ih2 ° 

corresponding to the coefficient of the first coordinate function In 
expanding the lefthand side of (10b), It can be shown that the expres- 
sion becomes Identically zero. On the other hand, the static equilib- 
rium of the composite plate requires that 



T2dz + 



(50) 


Thus, the first equation obtained from the series expansion of (10b) 
becomes an Identity, 0=0, and may therefore be Ignored. 

Now let us assume that the tractions are 

a^(z) = -a-j , (J2(z) = -02 » t^(z) = 0 , t2(z) = 0 , (51) 

where o-j and 02 are constant, and cos[iTk(z+h 2 )/(hi+h 2 )],(k= 0 ,l, 2 ,...) Is 
selected as the orthogonal system. By expanding OO) Into cosine series 
and considering the first N+1 terms we then obtain 

“n'^l^“n®^^^^''2°5n ^ln^®nk “n^Sn^’nk 
n*l 

* * °ln^Sk “n°5n*^nk * ^^'’l®5n * ^n^®nk 

“n^n^k ^ ^In^^nk ^ «n®5n^k^ = ° • 

(k=l,...,N), (52) 
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“"‘'’'"“in • ‘ “Sn^n^nlnk 

- * ''n‘>5r>"nk ' “n^BnVnk 


* (Vln - Zn*5n)"nk * \W„k 


■ ^’‘n'Sn * ^n®5n'5nk " “ti*5n*n*'nk’ 

“n *n k 

^ ^ [(■!) - cosY|^h2 + B(cosY|^h2-l )] 


= - a^nij^ - a2nj^ , (k=l N), 


(53) 


Re I a3{(x D, - Y„Cc„)a„^ + D. a X„b„„ 
n n In n on no on n n no 


n=l 


- (X C, + Y„D- )c„^ - a„Cc„X„d„„ 
n In n 5n no n 5n n no 


+ (X B, - Z A- )e + a Br X f 
' n In n 5n' no n 5n n no 


- (X A, + Z„B.„)g„„ - a A.„X„h„„} 

' n in n 5n ^no n 5n n no 

Do 

+ ^ (h] + eh2) ^ (h| - Bh|) = - - 02h2 


(54) 


where 


Xn = yV>*S7^l(V> • 

Y„ = (l+2v2)K^(a„a) + jLK,(o.^a) , 

Zn * (l+2vi)Ko(ona) + ^ Ki(ana) , (n=l,2...,N) , 

sinY^h- . 

("k * — , . Y|^ “ . (k=l,2..,N) , 
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1 COSYlHo 

1 sinY|,h, 1 cosYifhj 

♦ ^ s">(y|<%)I> 2 - ? V“t.>"2 - ’’ 

« slny^h^ 

-4 ^ sin(Yt-oi-)h, » (n»l,..,N; k»l,..>N) , 

i. Y|^-a„ K n c 

, COSYi,h« , COSYijhj 

»nk • - Vos(v«„)i'2 * TT7j;;^s'"(Yk%)'>2 


1 sinY|,h, 

? y’"(V“n)"2 


, sinYuh, 1 cosY|.h« 

, cosYth, 1 sinYt^hp 

- I (V|,-»nF ^'"^^k~°n>^2 * J V 7 -S 7 "2*'"<''k‘“n>''2 
, sinYuho 

- -K — ~ 2 [1 “ cos(Y|^"Otj^)h23 » (n*l,...»N; k*l,...»N) » 

(''k-“n> 

, cosYi,h, 1 slnY^hp 

'nk ° 1 T J ^ r V”n^^2 ^ I Y.+a„ ■ cos(Y|^%)h 2 ] 




, cosYfch, 1 sinYj^h^ 

®’"<Yk-»„)h2 ♦ C’-«‘>5(Yk-‘‘n>"23 • 

( n * 1,..,N; k * 1,..,N) , 

1 COSYi-ho 1 COSYuh- 

- J h2*'"< V“n>'’2 J Trj;^ [l-cos(r,,+a„)h 2 ] 

, sinYi-h, I sinYj^h- 

+ I VA » •’2“®(Yk+»„)''2 - ? (T|,t»Jk®'"<V“n>''2 


1 COSYth, 1 COSYt-hj , . , 

- ? h2*'"<’'k'“n>"2 * I ‘ “*<Yk-“n>''2l 


« sinY|.h^ \u ^ *iln( Yi )^o • 

♦ ^ -:^ hjCOsCYk-njhz - ^ 

(n*l»«»»Ni k*l».»>N) ♦ 

1 COSYi-hy / 4. \h 11 I 1 ^ y1"(ti*"..^N 

<■ COSYb^'o r / ^^l 11 I ^ Z. «nt )h- » 

* ’k n 

(k»l,2,..,N; n*l.....N)] . 

, coSYijhj j. \K + ^ sin(Yi,+o„)hi 

h,cos(Yk*a„)h, * J 5 ^''k n 1 

4 SiHYi,h2 • / 4. \h + i -r— tcOS(Yif*®n^^l”^ ^ 

^ ^k n 

Y \ t.,COk(Yk-a„)h, - \ Sin(Yk-«„)», 

. ^ h,sin(Y, ^ ’ 

(n=l»«»«N» k®l»-**»N) » 


COSYk^'2 • / + ^h * [cos(Yit'*'°n^^1~^^ 

sin(Y,^%)hl I Y,j^+a„ k n 1 


9nk*I Yk+«n 

, sin (Yk-a„)h, Y [COSlYkVr’l 

(n*1,...N; k"l,...H) • 


"nk* 


i'i!;^h,stn (Yk*»„)''i 

' ^k n 

1 / + \h - — T — s1n(Yt'*'°w)^'i 

^ \ ~ yj^ cos(Yk%)hi 2 lYk+On^ 
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n cosy^hp « cosykHa 

* 7 [cos(if|,-o„)h,-l] 

, sinv^h- , siny^h, 

* ? '’l “*<’'k‘“n>'’l - ? • 


k '‘n' 

(n«l,..,N; k«l,..,N) ; 


no 


‘'no 


"no 


no 


"no 


no 


^no 


no 


cosa^ho-1 
n c 


n 


hoCOSa^ho sino^h. 


& II & 

J. > » 

“n 

“n 

stna^hj 


a • 


n 


II 2 S 1 

(l-cosa^h2) 
^ ^ 

“n 

ijl 

n 

1-COSa^h^ 


“n 

> 

h,cos.„h, 

sina h, 
+ nl’ ' 

“n 

(X^ 

n 

sina^h, 


a_ * 


n 


ti,s1na„h, 

(l-cosoj^h^) 

n 



, (n»l,2..,N) . 


(55) 


and 8 is given by (42). Even though it is difficult to investigate the 
regularity of the algebraic system (52)-(54), the numerical results show 
very good convergence with increasing N. 


4. NUMERICAL RESULTS AND DISCUSSION 

In the numerical example considered the following material proper- 
ties and dimensions are used (Figure 1) 
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t 


* 0.3 


t 


^2 ■ 0.35 





( 56 ) 


The first 60 roots of the characteristic equation (43) obtained from the 
algebraic system (34)-(41) are given in Table 1. 

Two separate loading conditions are used to calculate the stresses. 
In the first it is assumed that 

~ "^2 

o ^/ o 2 * 1-964 , 0^21 


The second loading consists of a uniform pressure on the entire hole 
surface, namely 

(7pj,l(a»z) * 0^j.2(a»z) * -02 » 

^rzl^^*^^ * ~ ^ * (68) 

Tables 2-5 show the calculated results which are partially displayed 
also in Figures 2-5. Based on the calculated results one could make 
the following observations: 

(a) Away from the hole boundary generally the convergence is quite 
good. It becomes slower when the hole boundary is approached. For the 
loading given by (57) the discontinuity in traction a^.^. may be partially 
responsible for this. However, the main reason for the lack of conver- 
gence of the calculated results near the hole boundary appears to be 

”^he stress ratio 1.964 corresponds to (8a) and (8b) for the material 
pair under consideration. 
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the singular nature of the stress state at the intersection of the Inter- 
face and the boundary [1-3]. Thus* near the hole boundary more terms In 
the Infinite series were needed to obtain convergence comparable to that 
found In computing the stresses away from the boundary.* 

(b) As all stress components go to zero. However* the decay 
In Oggj Is much slower than that In and o^^j* (j»l*2). 

(c) From Figures 3 and 5 It may be seen that the relative magnl- 

tldues of the Interface stresses and are rather small. Also* 
the stresses corresponding to 0^*02 5) are an order of magnitude 

greater than those obtained for a-^ » 1.964 02 (Figure 3). A partial ex- 
planation for these results may be found If one considers the homogeneous 
plates separately. In a homogeneous plate axl symmetrically loaded by 
<j^^l(a*z) » -o-| * - Ot (0<z<hi) the stress state Is given by 

°rrl “ ■ (r7a)2 “ '’^eel ’ ®rzl * ° ’ °zzl * ° 


from which It follows that 
1 '^l 

‘esl ' 5U7 T777F ” '"rrl • 'rzl * ° ’ ' 22 I ° " 


(60) 


Thus* It Is seen that If the second plate Is axi symmetrically loaded by 

the two plates the displacements would be Identical along the Interface 
and the stresses and would be zero everywhere. In the example 
under consideration ui/yg® 2.077. Therefore* for a-j/og^ 1.964 one would 
expect the magnitude of the Interface stresses to be rather small. 
Similar observations may be made with regard to the comparison of 
and 0 ^^ In bonded and unbonded plates. On the other hand* for 0i*02 

*The numerical results given In the tables are obtained by using 20 to 
30 terms In the series for locations away from the hole and up to 60 
terms near the hold boundary. 
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one would expect higher Interface stresses because of the greater mismatch 
In the displacements along the Interface. 

(d) For the loading 0^/02 “1.964 since the solution Is close to that 

of a homogeneous plate, the thickness effect should not be significant. 
Indeed, by varying a/h-j It Is observed that the results do not change 
significantly. Also, In this case from Table 3 It may be seen that the 
variation of with 2 Is negligible, whereas for Table 5 shows 

a significant variation In a... Again, note that In the homogeneous 
plate Cgg Is Independent of z. 

(e) The calculated results Indicate that on the Interface 2>0 the 
stresses become unbounded as r approaches a, the hole radius. Theoret- 
ically, this Is Indeed known to be the case [1-3]. The solution given 
In this paper Is In terms of Infinite series, meaning that for z«0 and 
r»a certain series should be divergent. In problems such as this one 
would have to determine the eigenvalues in closed form for large 
values of n by examining the asymptotic behavior of the characteristic 
equation, determine the related eigenfunctions again In closed form, and 
try to separate and sum the divergent part of the Infinite series. Such 
a procedure seems to be quite Impossible for the problem under consider- 
ation. However, If one has a reasonably good solution for sufficiently 
small values of r-a, one may then follow an Indirect method to establish 
the singular behavior of the stresses in an approximate manner. To do 
this we note that from the plane strain solution of two bonded elastic 
quarter planes one may express the asymptotic behavior of the stresses 
for 2=0 and for small values of r-a as follows [3] 


’ 1 ] 




+ B 


ij 


l-a 

(r-a) + 


( 1 .j*r, 2 ) 


(61) 


where a Is the root of the characteristic equation In the strip 0<Re(o)<l 


(C0Sira + Cii + aCi2 Ci2C»(a+l )/ 2 ) (cosiro + C21 ’*‘^ 22 '* 

+ 0230(0+1 )/2) - (dii+d.|2o)(d2i+d22o) = 0 , (62) 
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where 


1 1*^1 ( 1 '*■'^'1 ) 3(1 -n^j) 

*=11 ‘I’ 2'(m‘+K") " * 

_ 6(1 -m^) 4(1 -m^) 

^12 " T+m^ » ^13 * “ l * 

_ 1 3(l-ni2) 

^21 * I ■ 2(iri^f) ■ * 


6(l-m2) 

^22 * l-hngKg * ^23 


4(1 -mg) 
l+IHgiCg 


3(1+Kg) 

^11 * J • ^(RnjKp • ‘*12 


1+IC^ 

l+mgtcg 


1+K^ 

mg+ic^ 


3 ( 1 + ic 2 ) l+Kg ^‘*‘‘^2 ^*'^2 

*^21 2Xnyh<^ ” 2(l-Hn^Ki ) * ^22 ” l-Hn-jK^ ' m-j+icg 


» Wg/*^! "*2 * *‘l^’^2 * "i “ ^-4v^ , (1»1,2) . 


For real material combinations it turns out that in 0<Re(o)<l (62) 
has only one root which is always real» and a»0 is not a root (meaning 
that there is no need to investigate the possible existence of a logarith 
mic singularity). For the material constants given by (56) and used in 
this paper a is found to be 

a - 0.048940 . (63) 

Thus, the approximate asymptotic behavior of the stresses around (z>0, 
r«a) may be established by assuming that a is known and by using the 
last two calculated points for in the expressions (61) to determine 
the corresponding constants A^j and The constants are usually 
referred to as the stress intensity factors. For the present problem 
they are found to be 
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A^j/lO'^Oga® 


1.964 

1.608 

0.776 

0.082 

-2.224 

1.000 

1.144 

0.562 

0.990 

1 

-20.918 

1 


This Is essentially a curve- fitting process to a smooth data. Conse- 
quently, for example, It was observed that the next point calculated 
from (61) Is rather In good agreement with the stresses given by series 
solution. 
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Table 1. The first 60 roots of the characteristic 
equation (43) 


n 

Re(o^h^) 


i"!) 

1 

2. 576640149. 

C 

1 

2 

1.293810176 


0.8298514402 

3 

2.885866996 

0.7551679918 


4 

4.384613298 


1.587138202 

5 

4.508247667 

0.9781794212 


6 

6.169990276 

1.228997641 


7 

7.t.82619850 


2.113409611 

8 

7.752063359 

1.340085172 


9 

9.342477506 

1.451922791 


10 

10.77125258 


2.453236163 

11 

10.92034586 

1.528341711 


12 

1 12.50012425 

1.602698086 


13 

13.94523374 


2.705188038 

14 

14.07574872 

1.661164351 


15 

15.65192294 

1.717452823 


16 

17.11030656 


2.905787940 

17 

17.22622651 

1.764831059 


18 

18.80069246 

1.810326620 


19 

20.26975009 


3.072577814 

20 

20.37415766 

1.850159284 


21 

21.94765603 

1.888423334 


22 

23.42541951 


3.215377334 

23 

23.52055621 

1.922785421 


24 

25.0934437 

1.955843539 


25 

26.57843058 


3.340252251 

26 

26.66594450 

1.986056608 


27 

28.23841688 

2.0515177619 


28 ; 

29.72949401 


3.451217584 

29 

29.81062486 | 

2.042135663 


30 

31.38280031 

2.068169865 


31 

3:. 87908457 


3.551071015 

32 

32.95478667 

2.092505940 


33 

34.52674116 

2.116052348 


34 

36.02753210 


3.641841058 

35 

36.09855559 

2.138231174 


36 

37.67034054 

2.159728360 


37 

39.17507330 


3.725045967 

36 

39.24201876 

2.180101156 


39 

40.81367039 

2.199880230 


40 

42.3218827 


3.801851958 

41 

42.38523874 

2.218718904 


42 

43.95678339 

2.23/036178 


43 

45.46809202 


3.873174424 
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Table 1 (Cont.) 


n 

te(a„h,) 

1 

44 

45.52826175 

2.254555535 


45 

47.09971910 

2.271613599 


46 

48.61380260 


3.939745157 

47 

48.67112277 

2.287986286 


48 

50.24250782 

2.303948118 


49 

51.7590938 


4.002158438 

50 

51.81384877 

2.319314977 


51 

53 38517321 

2.334313681 


52 

54.90402880 


4.060903496 

53 

54.95646096 

2.348791111 


54 

56.52773403 

2.362936830 


55 

58.04865835 


4.116387922 

56 

58.09897623 

2.376622132 


57 

59.67020534 

2.390007063 


58 

61.19302391 


4.16895489 

59 

61.24140822 

2.402982400 


60 

62.81259940 

2.415684503 
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Table 2 . Variation of the stresses at the Interface z«0, with 
r/a for 01/02 ■ 1.964 


r/a 




<J^2/(2a2Xl0‘^) 

1.025 

0.955 

0.460 

• 

-11.0 

1.037 

0.933 

0.449 

- 1.0 

-9.7 

1.050 

0.910 

0.438 

- 1.5 

-8.9 

1.075 

0.868 

0.417 

- 1.8 

-7.7 

1.100 

0.829 

0.399 

- 1.8 

-6.7 

1.150 

0.758 

0.365 

- 1.5 

-5.1 

1.200 

0.697 

0.335 

- 1.1 

-3.9 

1.300 

0.594 

0.286 

- 0.3 

-2.3 

1.500 

0.446 

0.215 

0.3 

-0.8 

1.750 

0.328 

0.158 

0.3 

-0.2 

2.0 

0.251 

0.121 

0.1 

-0.05 

3.0 

0.111 

0.054 

0 . 

0. 

4.0 

0.063 

0.030 

0 . 

0. 

5.0 

0.040 

0.019 

0 . 

0. 


Table 3. Distribution of stresses In 2 - direction for r/a*1.2 
and 0^/02 = 1.964 



2/h. 



°9el '^”2 


0. 

- 1.1 

- 3.9 

0.697 


0.25 

- 2.1 

- 2.1 

0.692 

1*1 

0,50 

- 1.5 

0.4 

0.687 

0.75 

- 0.5 

1.2 

0.682 


0.90 

- 0.1 

0.7 

0.679 


1.0 

0 . 

0 . 

0.677 


0. 

- 1.1 

- 3.9 

0.335 


-0.25 

0.1 

- 0.5 

0.339 

1*2 

-0.50 

- 0.4 

1.7 

0.344 

-0.75 

- 0.3 

1.6 

0.349 


-0.90 

-0.08 

0.9 

0.353 


-1.0 

0. 

0 . 

0.355 
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Table 4. Variation of the stresses at the Interface z«0, with 
r/a for 01/02 * 1.0 


r/a 


° 0 e 2 ^ 2‘’2 

02 ^/( 202 x 10 * 2 ) 

0 ^ 2 /( 202 ^ 10 * 2 ) 

1.025 

0.672 

0.324 

• 

- 10.4 

1.037 

0.653 

0.312 

- 0.9 

- 9.2 

1.050 

0.635 

0.301 

- 1.4 

- 8.5 

1.075 

0.605 

0.286 

- 1.7 

- 7.3 

1.100 

0.578 

0.273 

- 1.7 

- 6.4 

1.150 

0.529 

0.250 

- 1.4 

- 4.9 

1.200 

0.487 

0.231 

- 1.0 

- 3.7 

1.300 

0.415 

0.199 

- 0.3 

- 2.2 

1.500 

0.313 

0.151 

0.3 

- 0.8 

1.750 

0.230 

0.111 

0.3 

- 0.2 

2.0 

0.176 

0.085 

0.1 

- 0.05 

3.0 

0.078 

0.038 

0 . 

0 . 

4.0 


0.021 

0 . 

0 . 

5.0 

Hb9 

0.014 

0 . 

0 . 



Table 5. Distribution of stresses in z- direction for r/a*1.2 
and oi/oo * 1.0 


Ozzl/^^OgXlO"^) 



1=2 

0 . 

- 0.25 

- 0.50 


- 0.75 

- 0.90 

- 1.0 


- 1.0 

0.1 

- 0.4 

- 0.3 

- 0.07 

0 . 




® eei / 2 o 2 

- 

- 3.7 

0.487 

- 2.0 

0.444 

0.4 

0.399 

1.1 

0.350 

0.7 

0.321 

0 . 

0.302 

- 3.7 

0.231 

- 0.5 

0.270 

1.6 

0.318 

1.6 

0.369 

0.8 

0.400 

0 . 

0.423 
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Figure 2. Variation of the hoop stresses a 00 i(r,O) and ae 02 (r,O) with 
r/a for ai/az* 1-9M and z *0 





Figure Variation of the hoop stresses 0eei(r,O) and 0 ee 2 ('^»O) 
r/a for 01/03 * 1 



